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^ '■ ABSTRACT 

In the near infrared where detectors are limited by read-out noise, most inter- 
ferometers have been operated in wide band in order to benefit from larger photon 

^ \ rates. We analyze in this paper the biases caused by instrumental and turbulent 

effects to V 2 estimators for both narrow and wide band cases. Visibilities are 

^ | estimated from samples of the interferogram using two different estimators, V 2 

which is the classical sum of the squared modulus of Fourier components and a 
new estimator V 2 2 for which complex Fourier components are summed prior to 
taking the square. We present an approach for systematically evaluating the per- 
formance and limits of each estimator, and to optimizing observing parameters 
for each. We include the effects of spectral bandwidth, chromatic dispersion, scan 

. ^ | length, and differential piston. We also establish the expression of the Signal-to- 

Noise Ratio of the two estimators with respect to detector and photon noise. The 
V 2 estimator is insensitive to dispersion and is always more sensitive than the V? 
estimator. However, the latter allows to reach better accuracies when detection 
is differential piston noise limited. Biases and noise directly impact the dynamic 
range of reconstructed images. Very high dynamic ranges are required for direct 
exoplanet detection by interferometric techniques thus requiring estimators to be 
bias-free or biases to be accurately calibrated. We discuss which estimator and 
which conditions are optimum for astronomical applications especially when high 
accuracy visibilities are required. We show that there is no theoretical limit to 
measuring visibilities with accuracies as good as 10 -5 which is important in the 
prospect of detecting faint exoplanets with interferometers. 



Subject headings: Atmospheric effects — instrumentation: interferometers - 
methods: data analysis — techniques: high angular resolution — turbulence 
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1. 



Introduction 



Optical-infrared interferometers, are providing very high angular resolutions of order 
1 mas, a limit to be improved by the next generation instruments by a factor of 10. Several 
astrophysical topics topics require high dynamic range observations, one of which being the 
detection of exoplanets around their parent star. The required dynamic range scales with the 
ratio between the star and the planet, from a few 1000 for hot Jupiters to 10 10 for Earth-like 
planets in the near-infrared. The dynamic range in an image reconstructed from visibilities 
and phases can be approximated by the formula of Baldwin & Haniff (2002): 



where V is the visibility modulus and 5V the associated error, 5(f> is the phase error of the 
complex visibility and N vis the number of visibility points. Very large dynamic ranges are 
achieved with radio interferometers such as the VLA (120 000 for example in Hardcastle 
et al. 2003). However, infrared interferometers have a very limited number of telescopes 
with respect to radio and millimetric interferometers thus limiting the number of available 
visibilities. Eq. (1) shows that large dynamic ranges can only be achieved if visibilities are 
measured very accurately in this case. Other techniques to detect planets do not necessarily 
require images. Coude du Foresto et al. (1997) have shown that hot Jupiters could be 
detected by fitting a model to visibilities with absolute accuracies of 10 -3 . Future inter- 
ferometer facilities will have a sufficient number of telescopes to make better images. But 
for observations which require the highest dynamic range or for which high photometric 
precision is mandatory, model fitting will be the most rigourous technique. This process is 
sensitive to biases. At a more modest level of precision, possible biases must be considered 
even at the level of 10~ 2 when measuring visibilities which change substantially through the 
bandpass, e.g. spanning visibility nulls. 

Data reduction procedures must therefore focus on the use of estimators free from the 
well identified biases such as photon or detector noise. Optical-infrared visibilities are cor- 
rupted by biases induced by atmospheric turbulence. Atmospheric turbulence causes coher- 
ence losses which are difficult to calibrate because turbulence is not a stationary process 
and may vary between the source and the calibrator acquisitions. This important issue has 
been tackled by spatially filtering the beams using single-mode fiber optics, for example. 
Accurately calibrated data have been made available and allow the derivation of accurate 
astrophysical parameters of sources. Fibers, however, do not ameliorate temporal effects 
such as differential piston which remain a concern. 

As most interferometers have been operated in the near-infrared using small telescopes with 
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read-out noise limited detectors, most observations have been carried out in bands a few 
hundred nanometers wide to ensure sufficient signal. The current generation interferometers 
like VLTI or Keck benefit from large pupils but will still observe in wide bands at maxi- 
mum sensitivity. It is the purpose of this paper to study the sensitivity of wide band V 2 
estimators to parameters such as spectral resolution or instrument defects. Because piston 
is achromatic, the study of the sensitivity to turbulent optical path difference fluctuations 
applies equally well to the narrow and wide band cases. These biases must be taken into 
account when observations require maximum accuracy as is the case, for example, for direct 
exoplanet detection. Estimator bias must be taken into account in setting acquisition pa- 
rameters such as exposure time or fringe frequency and interferometer set-up such as spectral 
resolution. 

The general formalism used in this paper is presented in Section 2. The effect of spectral 
bandwidth on V 2 estimators is discussed in Section 3. We investigate the influences of longi- 
tudinal chromatic dispersion in Section 4 and of finite scan length in Section 5. We simulate 
the biases and errors caused by optical path difference turbulent fluctuations in Section 6. 

2. General formalism 

2.1. The interferometric equation 

We describe in this section the main equations required to produce the results presented 
in the next two sections. We have chosen the formalism of co-axial interferometers in which 
the optical path difference modulation is produced by moving optics (e.g. the Michelson 
interferometer) as opposed to multi-axial interferometers in which the optical path differ- 
ence varies with spatial coordinates in the image. The nature of the effects discussed in 
this paper is generally the same for the two types of interferometers. Multi-axial interfer- 
ometers are however more sensitive to optical path difference fluctuations and to bandwidth 
effects. The interferogram being multiplied by a non-flat window (the diffraction pattern for 
a two-aperture Young's experiment), fringe modulation degrades for non zero optical path 
differences when wide bandwidths are used. It is not the scope of this paper to discuss 
this issue and we will focus on co-axial interferometers for which the acquisition window is 
flat (the interferogram is multiplied by a temporal window of unit height). We choose a 
beamsplitter-type combiner with two complementary outputs. 

We note S(a) the spectrum of the source as a function of wavenumber a and F(a) the trans- 
mission of the detector filter or the spectrometer transmission for a given channel. Assuming 
perfectly balanced beams and perfectly corrected turbulence over each pupil, the intensity 
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Fig. 1. — Absolute difference between normalized and unnormalized V2 estimators for a 
100% visibility as a function of opd stroke. Calculation has been performed in the full K 
band for a G5 III star. Considering the current level of visibility accuracy of single-mode 
interferometers, the normalization factor can be dropped without causing a significant bias 
if strokes larger than 40 /im are used. 

of one of the two outputs can be written: 

r+00 

I(x = v.t)= / S(a)F(a)da (2) 
Jo 

r+00 

+ / S(a)F(a)V(a) 
Jo 

x cos (2%ax + 0(a) + 27rap(x)) da 

where x is the optical path difference (opd) between the two beams, t is time, v is the velocity 
of the optical path difference or fringe speed (twice the speed of a moving retro-reflector in 
a single-pass delay line system), V(a) is the visibility modulus as a function of wavelength, 
0(<r) is the phase of the visibility and p(x) is the differential piston. The first integral is 
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the total power collected by the interferometer. The second integral is the fringe pattern. 
Because of the finite bandwidth 5a of the filtering function F, the characteristic number of 
fringes in the interferogram derived from the coherence length equals a /5a. In the near- 
infrared with filter bands similar in width to the J, H or K photometric bands, the typical 
number of fringes is smaller than 10. 

In a perfect instrument, the total power can be accurately measured and the interfer- 
ogram can be normalized to provide a direct measurement of the real part of the mutual 
complex spatial degree of coherence 7 i2 (x): 

r+co 

7lf l (x) = / W(a)V(a) cos (2nax + <j>(a) + 2vap(x)) da (3) 
Jo 

where the weighting function W is: 

w(a) = jMm 

/ + « S(a)F(a)dc 

In the case of quasimonochromatic light (Goodman (1985)), YiT l ( x ) i s the real part 
of the product of the complex degree of coherence j(x) - which measures the temporal 
coherence of the source - by the complex coherence factor /i 12 which contains the spatial 
information on the source. From here on, the phase of visibilities is assumed to be zero (or 
7r when the sign of the visibility function changes) and 7 relates to the visibility through: 

r+00 

rif(x) = / W(a)V{a) cos (27iax + 27rap{x)) da (5) 
Jo 

/+00 
V e (a) exp (— 2ma (x + p(x))) da 
■00 

where we call V c (a) = \ (W(-a)V(-a) + W(a)V(a)) the extended visibility. When the 
piston term is equal to zero, the mutual degree of coherence 7 is therefore the Fourier 
transform of the extended visibility. Conversely, the Fourier transform of 7^ al (x) is the 
weighted visibility function versus wavelength, the principle of double Fourier interferometry 
(Mariotti & Ridgway (1988)). 



2.2. Visibility estimator 

2.2.1. Continuous signals 



In practice, the phase of interferograms is not accurately known hence it is not possible 
to measure the complex visibility. The observational goal is to estimate the modulus V of the 
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Fig. 2. — Difference between the monochromatic visibility modulus of a 43 mas uniform disk 
and the wide-band estimators as a function of spectral bandwidth at 2.2 /im. V\ is left, V-j 
is right. The peak errors in V\ occur at spatial frequencies corresponding to the visibility 
zeroes for the central wavelength of each filter bandpass. 



visibility. Let us first consider the piston- free case. In this special case, the interferometer is 
a generalized Fourier transform spectrometer working at a non zero baseline. It is possible 
to measure the weighted visibility modulus as a function of wavelength - the knowledge of 
the source spectrum then yields the visibility modulus as a function of wavelength: 



V{a) 



W{a) 



llT\ x ) ex P (2i7rax) dx 



(6) 



In presence of additive noise of variance S 2 an unbiased estimator is V 2 (a) = V 2 (a) — T? . 
This is why we only consider quadratic estimators in the following. 



2.2.2. Sampled signals 

Real instruments measure samples with spacing dx on a length Ac = N&c yielding the 
sampled mutual degree of coherence: 



7,W= ftx 



[ £)H 



/ X 

VAc 



(7) 



7i2 al ( x ) 1S t ne continuous real part of the complex mutual degree of coherence. The 

multiplication by the | | (f:) function expresses that it is sampled with a fix spacing between 

samples. The finite scan length is mathematically expressed by the multiplication by the 
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II (#0 function. This assumes that the time constant of the detector is negligible compared 
to the reciprocal of the sampling frequency. 

A discrete spectrum is obtained thanks to the convolution by the |_ J f^J function 
which periodically replicates the sampled temporal signal whose acquisition length Ac is 
finite. Replication is necessary to ensure that both 7^ and its spectrum V s are discrete. In 
addition, V s is also periodic. This is the classical mathematical expression of discrete Fourier 
transforms of sampled physical signals. 

For integrating detectors (most if not all infrared detectors) an extra convolution of 
7i2 al ( x ) by the door function II is required: 



We only consider the case for which sampling and integration times are equal. In 
visibility space, the spectrum is filtered by the sine function U(a&c) whose minimum value 
is 2/tc at Nyquist frequency. The impact of integration on the visibility estimator is to 
reduce the modulation and therefore to require a normalization by an ad-hoc factor. This 
factor depends upon the source spectrum and visibility distribution but also on the sampling 
frequency. The larger the sampling frequency with respect to the fringe frequency the smaller 
the contrast reduction and therefore the sensitivity of the estimator to spectral features. 
When fringe energy is spread over a large bandwidth because of vibrations or atmospheric 
piston, the filtering will remove energy and therefore reduce contrast. This effect appears 
when the source of bias is piston and is therefore considered in Section 6. Integration is 
otherwise not taken into account in the following for the sake of simplicity. 

Assuming that signals are sampled but not integrated, V s is equal to: 



7r 2 al (*) 



X 



X 



xn ( „ 

V Ax 



x 

Ax 



V S {<J) 



VJ(T)*n(aAx 



x 



|(o-Ar)l *|_ J(cr&) 



(9) 



where II, the Fourier transform of II, is a sine function of width 1/ Ax. The resolution 
of the measured visibility spectrum is therefore directly proportional to the interferogram 
length, a well-known fact in Fourier transform spectroscopy: R = a. Ax. In the following, 
the spectral information being replicated, we only consider the first N samples : 



V s (kda) 



V e (a)*Il(aAx) 



(kSa), k = 0, 



N 



(10) 



with da = 4-. 
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3. Effect of spectral bandwidth 



In this section we assume an infinite scan length and study the effect of spectral band- 
width which in practice cannot be negligible. The assumption of infinite scan length is equiv- 
alent to assuming that the source spectrum only contains spectral features whose width is 
larger than the Fourier spectral resolution Cases for which spectral features are smaller 
or close to the spectral resolution need to be modeled if high accuracy is required. 
This latter case is also a source of bias as illustrated by the following example. Suppose a re- 
solved source of continuum with a 50% flat visibility across the bandpass, and an unresolved 
spot in a 10 km/s line with the source and the spot emitting equal fluxes. The unresolved 
spot has a visibility of 100%. The line gets averaged with the continuum in the convolution 
by the spectral response of the interferometer. Assuming a scan length of 100 /im at 2 /im 
wavelength, hence a resolution of 50, the width of the line is 1/600 of a spectral element. 
The visibility being the linear average across a spectral element, it is of 0.50042 instead of 
0.5 with continuum emission only hence a bias of 0.08%. 

In the following we assume a smooth visibility spectrum and investigate the impact of band- 
width on estimators. 

With this assumption, Equation 10 has the simplified expression: 



We investigate the influence of the spectral bandwidth on the wide-band visibility estimators. 
Noise is neglected in the following but we use squared moduli to remain in the framework of 
an unbiased estimator in presence of additive noise. 

Two visibility modulus estimators can be built depending on whether real or imaginary parts 
of the extended visibility are summed prior to being squared: 



V s (kSa) = V e (kSa), k = 0, . . . , N - 1 



(11) 




iV-l 



(12) 



Eto 1 W(kSa)Sa 



nv-i 1 2 pv-i 



2" 



_fc=o J Lfc=o 



Both estimators have drawbacks and advantages. The first one, V±, is insensitive to 
atmospheric or instrumental differential phase effects such as chromatic phase dispersion 
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generated by air at visible wavelengths or by unmatched dispersion in single-mode fibers - 
it is for these reasons that it has been so widely used. Yet, it has to be normalized by the 
shape factor ^^Zq 1 W 2 {k5a)6a which depends on the source spectrum (Coude du Foresto et 
al. (1997)), hence possible bias effects must be considered. The second estimator, V 2 2 , is 
adapted from the multi-speckle mode visibility estimator of Berio et al. (1999). It is sensi- 
tive to differential phase and requires these effects to be negligible or at least stable enough 
to be calibrated. As V 2 2 applies a summation over linear quantities, the normalization factor 

r i 2 

~^2k=o W(k5o)5o is very close to 1 if long opd stroke scans are recorded. This makes this 
estimator much less sensitive to spectral features than V\. For example, the bias introduced 
by deep absorption bands is negligible in the case of Vi compared to V±. The normalization 
factor cannot however be neglected if short scan lengths are used as shown on Fig. 1. We 
have plotted the difference between V2 computed with and V2 computed without the normal- 
ization factor as a function of opd stroke. We have used the spectrum of a G5 III star in the 
K band. Considering the current level of visibility accuracy of single-mode interferometers 
of 0.2 % (Perrin et al. (2004)), the factor can be dropped without causing a significant bias 
if strokes larger than 40 /im are used in this case. The required opd stroke varies with the 
source spectrum. It is 60 /mi for a flat spectrum which is the most extreme example. When 
the filtered spectra of the calibrator and science target are similar, the normalization fac- 
tors disappear in the interfero metric calibration. Yet, the ratio of the normalization factors 
should be estimated when it is a potential critical source of bias. In the case of integrating 
detectors, the contrast loss due to the averaging of the fringe modulation on a finite time 
scale is directly taken into account by the normalization factor. For monochromatic signals 

it is directly equal to ^ ^[toIkj J with the particular value of 7r 2 /8 at four samples per fringe. 
The signal-to-noise ratios of both estimators are calculated in the Appendix. V\ is more 
sensitive than V2 by a factor 2 to 3 under realistic conditions. 



For both V\ and V2, the estimated visibility modulus is not exactly the object visibility 
modulus measured at the effective wavelength of the instrument. The reason is that the 
averaged squared visibility function in a band of finite bandwidth is different from the vis- 
ibility function. This is particularly true when the visibility function flips sign across the 
band. This effect has already been described in recently published papers for estimator V\ 
(Kervella et al. (2003); Perrin et al. (2004)). We have illustrated this with the example of 
a 43 mas uniform disk observed at a wavelength of 2.2 /im. We have plotted on Fig. 2 the 
difference between the monochromatic visibility modulus at the effective wavelength: 



V(S) 



2J 1 (tt0S) 

7T0S 



(13) 
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with S the spatial frequency and the uniform disk diameter, and the wide-band visibility 
modulus estimated with the Viand V2 estimators. Bandwidths range from 12.5 nm to 400 nm, 
the full K band. We have used flat spectrum and transmission. The difference between the 
two estimators is evident. V2 is very close to the monochromatic function even with the 
largest bandwidth. The difference is at most 0.6% and reaches maxima at the top of visibility 
lobes or at the inflexion point. The Vi estimator has the same features but the difference is 
much larger at the location of the nulls of the visibility function with a maximum difference 
of 4%. In both cases, the effect of bandwidth cannot be neglected for bands a few hundred 
nanometers wide. 




J I I L 



10 20 30 40 50 60 70 
Spatial Frequency (cycles / arcsec) 

Fig. 3. — Difference between the monochromatic visibility modulus of a 43 mas uniform disk 
and the V2 wide-band estimator as a function of dispersion strength in a 400 nm wide band 
around 2.2 /mi. The dispersion strength is expressed as peak-to-peak phase. 
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Fig. 4. — Bias of visibility estimators with respect to an infinitely long scan as a function of 
opd stroke. The source is a 43 mas uniform disk observed in a 400 nm wide band centered 
on 2.2 /im. V\ is left, V-j is right. 



4. Effect of chromatic dispersion 

When two interferometric beams are crossing different thicknesses of glass or air (at the 
blue end of the spectrum) longitudinal chromatic dispersion (or group delay dispersion) is 
produced. An additional phase term needs to be added to the fringe spectrum which takes 
the simple form: 

m = \D{o- o,f (14) 

where o"o is the center of the band and where the th and 1 st phase orders are set to (the 
interferogram is centered) . D is the dispersion coefficient (the second derivative of the phase) 
and can be expressed in units of /xradcm 2 . The reader is referred to Coude du Foresto et 
al. (1995) for more details on the notations and for a study of single-mode fiber dispersion 
for astronomical interferometry. The differential phase produced by dispersion has no effect 
on the V\ estimator by definition when scan length is infinite or when spectral resolution is 
large enough that phase is constant across a spectral element. In wide-band, the estimator 
gets biased by dispersion as is shown in Fig. 3 (we took a 400 nm bandwidth centered on 
2.2 /xm). We have chosen to represent cases for which the peak-to-peak phase produced by 
dispersion accross the spectrum is 1, 1.5 and 2 radians. Such cases can be encountered 
at short wavelengths because of the dispersion of the refractive index of air or when using 
fibers if arms are at unequal temperatures. A 1 K temperature difference over 1 km of fluoride 
fiber will produce a 10 mm difference of glass optical path length, with a 1.2 radian peak-to- 
peak phase at 2.2 /im (Kotani et al. (2004)). Any uncontrolled temperature variation will 
introduce visibility biases with the Vi estimator. Or equivalently at shorter wavelengths, 
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Fig. 5. — Simulation of the noise and bias produced by piston on a 43 and a 10 mas source 
observed in a 400 nm wide band centered on 2.2 fim. Two fringe frequencies have been 
simulated: 300 and 1000 Hz. See text for more details on piston sequences simulation. 



different air pathlengths will cause different dispersions which require a calibration for V^. A 
calibration can be achieved by measuring the Vi response on a point source or alternatively 
by measuring the dispersion law and correcting for it in the Vi estimator. However, the 
calibration is subject to changes hence to uncalibrated biases. 



5. Effect of finite scan length 



With temporal coding, data processing is very similar to that of Fourier transform 
spectrometers - after Fourier transforming the interferograms, the frequency samples have 
a resolution proportional to the scan length. As a consequence, the measured extended 
visibility spectrum is the convolution of the real extended visibility by the Fourier transform 
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spectrometer response as expressed in Equation (10). Because of the convolution, measured 
Fourier components are the weighted sums of spectral elements. The V±, for example, is 
therefore not the pure quadratic sum of spectral elements. A bias is therefore to be expected 
at spatial frequencies where the visibility varies quickly or flips sign. This is illustrated in 
Fig. 4 where the biases of the visibility estimators for a 43 mas uniform disk observed in a 
400 nm wide band centered on 2.2 /im with respect to infinitely long scan estimators have been 
plotted for different opd stroke values ranging from 256 to 2 /mi. When the stroke is 2 /im the 
estimator at 2 /im wavelength is the famous ABCD estimator. As expected, the V\ estimator 
bias is large where the visibility function of a 43 mas uniform disk comes to a null. The bias 
becomes lower to the best current accuracy of 0.2% for scan lengths longer than 64 /im. The 
V2 estimator naturally turns out less sensitive to the effect as the spectral components are 
first linearly co-added before squares of the imaginary and real parts are summed. Except 
for very short scans, the bias is always negligible. This calculation shows however that the 
estimator needs to be properly designed to not introduce a bias in the interpretation of 
visibilities if short scan lengths are preferred to longer ones. This particularly holds for the 
multi-axial combination for which a few fringes across the diffraction pattern are measured. 



6. Effect of differential piston 

Unlike dispersion effects or purely acquisition or instrument dependent biases of the 
estimators, errors due to piston cannot be calibrated so easily. Piston errors have coherence 
times of a few milliseconds and their characteristics are not stationary. The best solution 
is to include a fringe tracking system that will correct for the lower frequency part of the 
piston. Still, with such a system, only piston errors in the band pass of the fringe tracker 
are suppressed. Higher frequency errors will blur the fringes during exposures and reduce 
the contrast - a source of bias if the fringe tracker performance on the source and on the 
calibrator are different or if turbulence has evolved between measurements. An alternative 
is to scan fringes fast enough that piston errors can be frozen during an acquisition at the 
cost of sensitivity, which is equivalent to benefitting from a fringe tracker with a high band 
pass. We investigate here the influence of piston in this latter case. 

We first have chosen to simulate piston and data acquisition sequences which corre- 
spond to real data acquired on the supergiant a Orionis and its calibrators and reported in 
Perrin et al. (2004). We have therefore considered two diameters: 43 and 10 mas and two 
fringe frequencies (v/X with the notations of Section 2): 300 and 1000 Hz. Typical seeing 
conditions are on the order of 1 arcsec with a wind velocity of lOm/s giving ro = 45 cm (i.e. 
exactely the size of the 45 cm telescopes) and a coherence time for turbulent phase of 45 ms 
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in the K band. We have simulated scan lengths of 64 /im with 4 samples per fringe. For the 
fringe frequencies of 300 and 1000 Hz, 99 and 29.7 /im of opd are respectively scanned in a 
coherence time of 45 ms. The piston is therefore mostly frozen with the higher frequency and 
only partially with the smaller one. However, in the latter case, the whole 64 /im acquisition 
length is required to avoid biases due to finite scan length. Piston stability is primarily re- 
quired over the coherence length which contains most of the fringe power. No additive noise 
has been injected in the simulations - only pure piston noise is considered here. Visibilities 
have been measured for each simulated interferogram using the V\ and Vi estimators. 1000 
interferograms are processed for each spatial frequency sample. The piston sequences are 
computed using the method of Perrin (1997). The baseline used to compute the piston rms 
and power spectral density varies with spatial frequency according to: B = X.S. The lower 
the spatial frequency the smaller the piston rms. This is a little pessimistic for the longest 
baselines as recent studies have shown that the outer scale of turbulence at current best sites 
is about 20 m (see e.g. Martin et al. (2000)) inducing a saturation of the rms piston to 
less than 20 /im (Davis et al. (1995)) or equivalently at baselines of order 20 m depending 
on models. In our simulation the piston rms does not saturate and linearly increases with 
baseline up to 35 m. Note that the rms of piston is at spatial frequency. 

The result of the simulation is presented in Fig. 5. The differences between unpistoned 
and pistoned estimated visibilities are plotted. The error bars are the standard deviations 
computed from the 1000 interferograms. Estimates are biased when differences are more than 
3 a from 0. We conclude from the simulation that V\ is far less biased by piston than the V<i 
estimator. At the level of our simulation, piston can be considered a zero-mean noise for the 
Vi estimator, biases being smaller than 0.1% at 300 Hz fringe frequency. This is intuitively 
reasonable as the V\ estimator is not sensitive to interferogram shifts and by extension poorly 
sensitive to the first orders of piston whereas the V2 estimator being a coherent estimator is 
sensitive to all piston orders including the lowest ones. 

Piston is clearly a bias for V2 at low fringe frequency. Interestingly, the statistical scatter 
in the simulation for V2 is always much smaller than for the V\ estimator - perhaps thanks 
to the coherent averaging of spectral elements. 

At higher fringe frequency however, when piston is almost completely frozen, the bias 
on V2 is usually far less or less than 0.01%. In practice, only bright sources in good weather 
conditions can be observed at high fringe frequency. With the parameters used for our 
simulation, the statistical error on visibilities generated by piston is of 0.1% or larger if less 
than 1000 scans are acquired. This is in accord with current limits in visibility accuracies 
measured by single-mode interferometers which can thus be considered unbiased by piston 
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Fig. 6. — Simulation of the noise and bias produced by piston on estimators for a point-like 
source. See text for parameters. The insets are close-up views of the curves between 500 
and 2000 Hz. Integration is taken into account for the right panel whereas signals are just 
sampled in the left panel. 

if suitable fringe frequencies are chosen - as in, for example, Perrin et al. (2004). 

It turns out from this discussion that the two estimators have different properties de- 
pending on fringe acquisition parameters. We have computed the bias and the statistical 
error for the two estimators as a function of fringe frequency for the following set of pa- 
rameters (Fig. 6, left panel): 1 arcsec seeing, infinite outer scale, 10 m/s wind velocity, K 
band, 64 fim scan length, 4 samples per fringe, 45 cm telescopes, unresolved source (100% 
visibility), 40 cycles/arcsec spatial frequency (18 m baseline), 10000 interferograms. The 
following conclusions do not depend upon bandwidth and also apply to narrower bands. 
The V\ estimator is biased at more than 1% below 100 Hz fringe frequency. Same for Vi 
but below 300 Hz. V\ is unbiased above 300 Hz. Above 800 Hz, V2 becomes much more 
accurate than Vi and reaches accuracy levels below 10~ 5 . As a conclusion, V\ is better suited 
for lower fringe frequencies than V2 but V2 provides much better accuracies at larger fringe 
frequencies or with a fringe tracker with a band pass larger than 800 Hz. 

The estimators have also been simulated taking integration into account between con- 
secutive samples. The general tendency is the same except at frequencies below 200 Hz for 
which the bias is more important. The integration slightly blurs the fringes if the fringe 
frequency is too low. The effect is more important for than for V 2 2 which performs a 
linear sum of Fourier components. Results are displayed in the right panel of Fig. 6. 

These conclusions can be extended to other parameter values. First, the noise and bias 
produced by piston are multiplicative. The piston effect can be scaled from Fig. 6 for any 
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visibility. The fringe frequency axis can be rescaled as a function of wind speed according 
to: / x ( 10m / s ) , We have chosen a baseline of 18 m which is supposed to be just below the 
saturation of the piston rms. Longer baselines should therefore generate the same amount 
of piston whereas shorter baselines will have less piston making this simulation a little bit 
pessimistic. Scaling this graph to better or poorer seeing is unfortunately tricky as the link 
between visibility errors due to piston and seeing is not linear. In first approximation this 
graph can be considered for any seeing value around 1 arcsec. The size of telescopes and 
the baseline will change the cut-off frequencies of the piston spectrum whose consequence is 
difficult to guess from Fig. 6. In the same vein as for seeing, results from this graph should 
be considered applicable in first approximation. New simulations should be considered for 
significantly different values for these parameters. 

This study shows that the fringe frequency should therefore be tuned with respect to 
the piston conditions and the accuracy to be achieved in the same way as the number of 
exposures is set in order to achieve the required signal-to-noise ratio. The visibility estimator 
should also be chosen accordingly. 



7. Conclusion 



We have shown that visibility estimators must be computed with care for low spectral 
resolution measurements. We have compared the performances of two visibility estimators 
under varied instrumental or turbulent conditions. The V\ estimator tends to be more 
robust to all effects. It is even intrinsically insensitive to longitudinal dispersion. The Vi 
estimator can be used if differential piston is frozen and if dispersion effects are either stable 
or measurable with great accuracy. Besides, the spectral slope of the phase has to be removed 
accurately. Vi is more sensitive to additive noise than V\ by a factor 2 to 3. A first important 
conclusion is that the biases of the V\ estimator are smaller than the best accuracies (0.2%) 
measured for fringe frequencies larger than 200 Hz. A second important conclusion is that, 
although the V2 estimator is more sensitive to phase defects, it has the potential to reach 
accuracies as good as 10 -5 for high fringe frequencies or with a fringe tracker with a large 
bandpass. If the visibility accuracy is limited by piston rather than additive noise then the V2 
estimator is better. When the observation is performed with the goal of achieving excellent 
accuracy, effects should be simulated to tune the acquisition parameters such as the length 
of the interferograms, the acquisition frequency and the number of recorded exposures. High 
dynamic range detections are therefore achievable at this price and with suitably designed 
instruments. 
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A. Estimator signal-to-noise ratio 

The signal-to-noise ratio (SNR) of the ABCD estimator has been derived by Tango & 
Twiss (1980) for the photon noise case and by Colavita (1999) for both photon and detector 
noise. We present here a derivation of the SNR for the V 2 and V 2 estimators for any number 
of samples. 



A.l. Sampled interferogram formalism 

Notations are derived from the general interferogram expression of Eq. (2): 

L — Ci + Mi i = . . . JV - 1 



(Al) 



in which A" is the number of samples, Cj and Mi are respectively the continuous (unmod- 
ulated) and modulated parts of the interferogram in photon units. These quantities are 
corrupted by both source photon and detector noise - the realization of the latter is noted 
bi for each i. Noise-free quantities are noted (ensemble average): 



Ci 

mi 



(Mi) 



(A2) 



In the following, we assume that the subtraction of the continuous part of the interferogram 
and its normalization to get the expression of estimators as in Eq (12) are noise- free processes. 
The justification is that the estimate of the continuous part can be low-pass filtered and 
contains more photons than the modulated part hence a larger signal-to-noise ratio. Let 
Bk, k = ... A" — 1 be the discrete spectrum of the modulated samples: 



Bi 



N-l 

E 

i=0 



-2iirik/N 



(A3) 



The unbiased estimators are proportional to the following expressions (taking into account 
proportionality factors does not change the final expression of the SNR): 



N-l 



N-l 



V 2 oc 



V4 



E i^i 2 - 

;=0 
N-l 



N 2 ~ 2 



a 



k=0 



^E 

i=0 



(A4) 



N-l 



k=0 



N 2 a 2 



N 2 m 
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In the following, in order to make the reading easier, the oc signs are replaced by equal signs. 
The expectations of the estimators are: 

N-l 



i=0 



(V 2 2 ) = N 2 m 2 

For the V% estimator not to be biased it is therefore very important that the th sample be 
centered on the white light fringe or, equivalently, the slope of the spectral phase must be 
accurately removed. 



A. 2. Variance of V 2 

The bias terms are dropped as they do not contribute to the variance. As a consequence: 

JV-l N-l 

var^ 2 ) = ^var(| J B fc | 2 )+ ^ cova^f , l^ 2 ) (A6) 

k=0 k^q;k,q=0 

which can be rewritten as: 

var(Vf) = ^covardfifcl 2 , \B q \ 2 ) (A7) 

k,q 

All indices run between and N — 1. Working on individual terms we shall first derive: 

covard^l 2 , \B q \ 2 ) = (\B k \ 2 \B q \ 2 ) - {\B k \ 2 ) {\B q \ 2 ) (A8) 



The expansion of the product \B k \ \B q \ yields: 

\B k \ 2 \B q \ 2 = M j M i M rnM n e- 2tn{ij - l)k/N+{m - n)q/N) (A9) 

j,l,m,n 

The second product in the covariance can be expanded as well providing the expression: 
covar(|5 fc | 2 , \B q \ 2 ) = ^ ((MjMiM m M n ) - (MjMi) (M m M n ))e~ 2 ^- l ^ N+ ^- n ^ 

j,l,m,n 

(A10) 

and eventually: 

varfX 2 ) = ^covard^l 2 , \B q \ 2 ) (All) 

k.q 

= ((MjMiMmMn) - (MjMi) (M m M n ))e- 2m{ij - l)k/N+im - n)q/N) 

k,q j,l,m,n 
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The j, l,m,n and k, q indices being independent, sums can be switched in whatever order in 
the above expression: 

var (V r 1 2 ) = J2 J2(( M i MlMmMn > ~ ( M i M i) (M m M n ))e- 2znij - l)k/N ^ e -*«{rn-n) q /N 

j,l,m,n k q 

(A12) 

The sum over q equals zero except when m = n in which case it equals N. Thus : 

vax(V?) = NJ2 E« M ^0 " < M i M <> «» J2 e ~ 2l7rU ~ l)k/N ( A13 ) 

j,l,m k k 

In the same way, the last sum equals zero except for j = I and the expression of the variance 
simplifies to: 

var(V?) = N 2 J2 (M]Ml) - <M 2 ) (M 2 m ) (A14) 

Individual terms are equal to zero when j ^ m as additive noises are uncorrelated yielding 
a still simpler expression: 

JV-l 

var(^) = A 2 £ {Mf) - <M 2 ) 2 (A15) 

i=0 



A. 3. Variance of V 2 2 

The case of is much simpler as one may easily show that the expression of the 
unbiased estimator is: 

V 2 = N 2 M 2 - N 2 a 2 - N 2 m (A16) 

As a consequence: 

var(V 2 2 ) = A 4 var(M 2 ) (A17) 



A. 4. Detector noise 

Modulated samples are considered only and are written: 

Mi=mi + b, (A18) 

The variance of read-out noise is a 2 . The fourth and second moments of the gaussian 
distribution Mj are easily calculated: 

(Mf) = m\ + 6m 2 a 2 + 3a 4 (A19) 
<M 2 > = m 2 + a 2 
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and the variances of the estimators are: 

var(Vf) = 2iVV 



N-l 



2 m l + N ^ 



i=0 



(A20) 



var(V 2 2 ) = 2iVV [2m 2 + a 2 ] 



A. 5. Photon noise 

Both the continuous and modulated signals contribute to photon noise. The number of 
photons to be considered to calculate the Poisson fluctuation per sample is therefore c$ + mi 
which we note rij. As in Tango & Twiss (1980), we consider that the statistics of different 
samples for broadband thermal radiation are uncorrelated. Consequently, the fourth and 
second moments of the Poisson distribution M,- are: 



{Mf) 



Hence the variance of the estimators: 



var(V 1 2 



n 2 + ni 



N-l 
i=0 



(A21) 



(A22) 



var(V 2 2 



N 4 n (4n 2 + 6n + 1) 



A. 6. Variance of estimators with filtered samples 

In practice, the fringe peak covers only a limited range of the Fourier components 
which are summed to compute the estimated visibility. In the derivation above, all Fourier 
components have been used thus increasing the variance. 

Let us assume that the fringe peak has a width of Ato which corresponds to the spectral 
bandwidth. Whatever the nature of the noise, source photon or detector noise, its Fourier 
transform is a white noise and therefore the squared modulus of its Fourier transform is 
still a white noise. Selecting the spectral window of width Ato for the computation of the 
estimator amounts to reducing the variance of the noise by a factor At/At , with At = l/&c. 
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A. 7. Signal-to-noise ratios 



Signal-to-noise ratios can now be derived: 
SNR X = 



At 



S^N-l 2 



SNRo 



At 



At 



2a 2 



(A23) 



JV-l 



Ar ^a 2 [2m 2 + o- 2 ] + n (4ng + 6n + 1) 



It is interesting to see that the SNR of the second estimator is independent of the 
total number of samples. The SNR of the first decreases with the number of samples. In 
practice, unless N tends towards infinity, V 2 always has a larger SNR than ■ Although 
it has been derived differently, the V 2 SNR estimate yields comparable results to the SNRs 
published by Tango & Twiss (1980) and Colavita (1999). An extra term however needs to 
be added to the result of Colavita (1999) and the variance due to detector noise should read 
2 (N) V 2 a 2 + 16a 4 instead of 16a 4 only ((N) is the average total number of photons and 
fringe phase is constant during a sample integration). The derivation of SNRs is illustrated 
by two examples in the detector noise and photon noise limited regimes. 



Detector noise limited regime Let us consider an interferogram with 128 samples in a 
400 nm wide band centered on 2.2 /im with a sampling &c = 0.5 fim. The filter has a uniform 
response across the band and the spectrum is flat. We assume the amplitude of the white 
light fringe is mo = nV with n a number of photons. We choose V = 1. The detector noise 
standard deviation is set to a = n x 0.2. With these particular values, the two SNRs are: 

SNR 1 = x 7.63 (A24) 

AaT 

SNR 2 = J— x 2.48 
V At 

The V 2 estimator is therefore better than V 2 for these realistic conditions which is the 
case as long as the number of samples remains smaller than a few thousands at high noise. 



Photon noise limited regime The same parameters are used for the photon noise limited 
regime. The number of uncorrelated photons per sample is set to n = 10 yielding the SNRs: 
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, Ax 

SNRi = x j— x 1.33 (A25) 



SNR-2 = \l— x 0.54 

ZVJn 



The same conclusion applies as for the detector noise limited case. is better than V^ 2 
as long as the number of samples remains smaller than a few thousands which in practice is 
not required. 



